Disordered Kondo Nanoclusters: Effect of Energy Spacing 
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Exact diagonalization results for Kondo nanoclusters alloyed with mixed valence impurities show 
that tuning the energy spacing, A, drives the system from the Kondo to the RKKY regime. The 
interplay of A and disorder gives rise to a A versus concentration T = phase diagram very rich 
in structure, where regions with prevailing Kondo or RKKY correlations alternate with domains of 
ferromagnetic order. The local Kondo temperatures, Tk, and RKKY interactions depend strongly 
on the local environment and are overall enhanced by disorder, in contrast to the hypothesis of 
"Kondo disorder" single-impurity models. 

PACS numbers: 75.20.Hr, 75.75.+a,75.40.Cx,75.40.Mg 



Magnetic impurities in non magnetic hosts have been 
one of the central subjects in the physics of strongly cor- 
related systems for the past four decades^]. Such endur- 
ing, ongoing research effort is motivated by a constant 
shift and increase of scientific interest over the years, from 
dilute to concentrated impurities Q , from periodic Q 
to disordered samples ii, and from macroscopic '?! to 
nanoscale phenomena [8| . 

For macroscopic strongly correlated d- or /- systems, 
the central problem is understanding the variety of mag- 
netic phases observed experimentally upon varying the 
pressure and species concentration, often giving rise to 
non-Fermi-liquid (NFL) behavior at low temperatures 
[1, 0, 0, 0- Various disorder-driven models have been 
proposed to explain the experimentally observed^ NFL 
behavior|^ 1^ 0, 0- The phenomenological "Kondo dis- 
order" approaches Hoj , based on single-impurity mod- 
els, assume a distribution P{Tk) of Kondo temperatures 
Tk = exp(-l/ p(Ef)J, caused by a distribution of ei- 
ther /— electron-c-electron hybridization or of impurity 
energy levels. These models rely on the presence of cer- 
tain sites with very low Tk spins leading to a NFL behav- 
ior at low T. An open issue in such single-site methods is 
whether the inclusion of RKKY interactions would renor- 
malize and eliminate the Iow-Tk spins. An alternative 
view is the formation of large but finite magnetic clusters 
(Griffith phases) within the disordered phase through the 
cornpetition between the RKKY and Kondo interactions 
0, E|. The first interaction being responsible for the 
quenching of the local /-moment (LM) via the screen- 
ing of the conduction electrons, whereas the latter being 
responsible for magnetic ordering. 

On the other hand, the relevance of small strongly cor- 
related systems to quantum computation requires under- 
standing how the infinite-size properties become modi- 
fied at the nanoscale, due to the finite energy spacing 
A in the conduction band H E [H H [H. Recent 
advances in nanotechnology have made possible experi- 
ments in extremely small samples, stimulating a resur- 
gence of interest in the single-impurity Kondo physics at 
the nanoscale. For such small systems, controlling Tk 
upon varying A is acquiring increasing importance since 



it allows to tune the cluster magnetic behavior and to 
encode quantum information. While the effect of size 
or energy spacing on the single-impurity Anderson or 
Kondo model has received considerable attention recently 
H> E> El E> El , their role on dense strongly correlated 
clusters with or without disorder remains an unexplored 
area thus far. 

In this work we present exact diagonalization calcula- 
tions for d- or /-electron nanoclusters to study the ef- 
fect of disorder and energy spacing on the interplay be- 
tween the Kondo and RKKY interactions. While it is 
well known that the cluster properties depend on cluster 
geometry and size|l6|. the motivation of the present cal- 
culations is that they treat exactly the Kondo and RKKY 
interactions and they provide a distribution of local Tr-'s 
as renormalized by the presence of intersite /-/ interac- 
tions. Our results show that i) tuning A can drive the 
nanocluster from the Kondo to the RKKY regime, i.e. a 
Doniach phase diagram in small clusters; ii) the A 
versus alloy concentration T — phase diagram exhibits 
regions with prevailing Kondo or RKKY correlations al- 
ternating with domains of ferromagnetic (FM) order; and 
iii) the local Tr-'s and the nearest-neighbor RKKY inter- 
actions depend strongly on the local environment and 
are overall enhanced by disorder. This disorder-induced 
enhancement of Tk in the clusters is in contrast to the 
hypothesis of "Kondo disorder" models for extended sys- 
tems. 

We consider a random binary alloy cluster, An^^Bx, 
of N=6 sites and different number of B atoms, x = Q- 
N, arranged in a ring described by the half-filled (Nei — 
12) two-band lattice Anderson Hamiltonian; in standard 
notation, 

ij G i<T i 

A. Effect of disorder 

We introduce binary disorder in the /-orbital energy 
(ej or ) and in the intra-atomic Coulomb energy U 
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Ua or Ub), to model two different types of atoms: a 
Kondo-type A atom with = —U^/2= -3 (symmetric 
case) and a mixed- valent (MV) type B atom with e_B = 
-2 and = 1. Both types of atoms have the same on- 
site hybridization Va — Vb — 0.25. For t = 1, this choice 
of parameters leads to a degeneracy between the doubly- 
degenerate c-energy levels, ek — ~t and the energy level 
eJ + Uoi the MV atom. Upon filling the single particle 
energy levels for any x, N — x (x) electrons fill the 
{ej) levels , and two electrons the -2t conduction energy 
level, with the remaining N — 2 electrons accommodated 
in the a; -1-4 degenerate states at ~t. This in turn results 
in strong charge fluctuations. 



The configurations for a; < 3 are shown in Fig. f, 
left panel, along with the value of the spin, Sg, of the 
ground-state. The A (B) atoms are denoted by closed 
(open) circles, respectively. Except for the homogenous 
cases (a;=0 and x = 6), with a Sg = ground state, for 
all X there are configurations with Sg ^ 0. The aver- 
age occupation and average LM for the periodic Kondo 



and MV lattices are < >= 



1, < 
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0.99, 



and < > = 1.6, < (M/ ) > = 0.43, respectively. 
We carry out a detailed analysis for x=l (Sg =2) to 
demonstrate the FM transition induced by a single MV 
atom in an otherwise Kondo cluster. Studies of extended 
systems have reported similar occurrence of ferromag- 
netism in the MV phase^j^l- As expected, the singlet 
ground state of the x = Kondo cluster is character- 
ized by n.n. anti- ferromagnetic (AF) /-/spin correlations 
(< Sf{i)Sf{i + !)>=- 0.58). The introduction of 
a MV atom renders them ferromagnetic. Since Ub is 
small, the B impurity tends to remove charge from the 
the conduction band, in particular from the A:-state with 
Cfe = — i, which has large amplitude at the B site and at 
the opposite A site across the ring. Such a depletion is 
different for the two spin states, thus yielding a maximum 
value for the f-moment of the MV atom. The /-/ spin 
correlation function between the Kondo and MV atoms 
are AF (< Sf{i)Sf{i + !)>=- 0.23), while they are 
FM among the Kondo atoms (< Sf{i)Sf{i -|- 1) > = 
-1-0.94). A similar result was recently found in ab ini- 
tio calculations 20], where introducing a nitrogen impu- 
rity in small (1-5 atoms) Mn clusters induces ferromag- 
netism via AF coupling between the N to the Mn atoms, 
whilst Mn-Mn couple ferromagnetically. We find that 
there is a crossover in Sg from — > 1 ^ 2 — > (Fig. 1, 
right panel) indicating a reentrant nonmagnetic transi- 
tion around = 2. This almost saturated FM Sg = 2 
domain is robust against small changes in Ub, V, ^a, Ua, 
cluster size {N = 7), and band filling (A^e; — 10) provided 
that the Kondo atom has a large LM. 



We have also calculated the effect of disorder on the 
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FIG. 1: Left panel: Alloy configurations for various concen- 
trations a; < 3 (the a; > 3 cases are obtained by exchanging 
closed and open circles). For each a; < 3 configuration, the 
value of the ground-state spin Sg is reported. Right panel: 
Energy difference (in units of 10~^t ) between the lowest S < 2 
eigenstates and the ground state as function of es. 
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FIG. 2: Temperature dependence of the average f- 
susceptibility for different alloy concentrations. The inset 
shows the low-temperature behavior. 



temperature-dependent average susceptibility, xlC^), 
kBTxliT) 



< Sf{{)S^"' >«^^ .(2) 



Here, 5*-^°* is the z-projection of the total spin (both 
/- and c-contributions) , ac^ are the configuration- 
dependent exact many-body states, and J2ca 
denotes exact configuration averages. In Fig. 2 we 
present TxliT) as a function of temperature for differ- 
ent X. As T -> (inset Fig. 2) Txl{T) approaches a 
finite value for x — I — A while it vanishes exponentially 
for x—0, 5 and 6. This is due to the fact that the 
former concentrations involve some configurations which 
are magnetic, while the latter have singlet ground 
states (Fig. 1). The stronger (weaker) low-temperature 
dependence for x = 1 {x — 2 — 4) is due to the smaller 
(larger) spin gap between the ground state and the 
lowest excited states. The magnetic susceptibility 
displays also a magnetic crossover upon varying x, and 
reveals a Curie-like divergence at low T for a; = 1 — 4. 
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TABLE I: Local Kondo f-c and n.n. /-/ spin correlations 
functions and the local Kondo temperature (in K) for two 
values of t (in eV). The average energy spacing is A = 4t/{N — 
1) = 4f/5. 



a) t=0.05 



b) t=1.0 





< Sf{i)S^ii) > 


< Sf{i)Sf{i + l) > 


TK{^) 


t=0.05 


-0.626 


-0.322 


551.8 


t=1.00 


-0.011 


-0.584 


173.4 



The temperature-dependent results for the specific heat, 
not reported here, show corroborative evidence of this 
disorder-induced magnetic crossover. 

B. RKKY versus Kondo: Effect of energy spacing 
Next we address a number of important open issues, 
namely (1) the effect of A on the interplay between 
RKKY and Kondo interactions in disordered clusters, (2) 
the characterization of the single- impurity "Kondo cor- 
relation energy" Tr- in a dense- impurity cluster and (3) 
the effect of disorder and A on the distribution of the 
local Tfj-'s. In the following, es — —2. 

In contrast with previous studies, which introduced a 
phenomenological distribution P{Tk) of single-impurity 
Kondo temperatures, the advantage of the present cal- 
culations is that one calculates exactly the Kondo cor- 
relation energy: we employ the so-called "hybridization" 
approach |2lj. with Tr- defined as 



ksTKil) = Eg{V, = 0) - Eg, 



(3) 



where Eg{Vi — 0) is the ground-state energy of the dense- 
impurity cluster when V is set to zero at the ith site. 



Eq.(3) reduces to keTK = Eband-EF + ef -Eg\]J,\2^ in 
the single impurity case. Here, Ep is the highest occupied 
energy level in the conduction band and E^and is the 
conduction band energy. This definition of the local Tk 
takes into account the interaction of the /-moment at 
site i with the other /-moments in the system |23| . 

In Table I we list for the periodic, a;=0, case the local 
Kondo f-c spin correlation function < S^{i)S^{i) >, the 
n.n. /-/spin correlation function < Sf{i)Sj-{i -\- 1) >, 
and the local Kondo temperature for two different values 
of t (The energy spacing is A = 4i/{N - 1) = 4i/5). 
As t or A decreases the f-c spin correlation function is 
dramatically enhanced while the /-/ correlation function 
becomes weaker, indicating a transition from the RKKY 
to the Kondo regime. This is also corroborated by the 
increase in the local T^(i). The energy spacing affects 
not only the magnetic (A) atoms but the MV atoms as 
well. Thus, increasing t drives the B atoms from the 
non-magnetic, NM [nj k, 2), to the MV and finally to 
the Kondo regime. 

In Fig. 3 we present the local TV (i) as a function of the 
local f-c spin correlation function < Sf{i)S^{i) > for all 
Kondo (A) atoms in the singlet ground state at any con- 
centration X for t= 0.05 and 1.0. Note the different scales 
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FIG. 3: A-atoms: Local Kondo Temperatures (in K) vs the 
local f-c spin correlation function, for different configurations 
and two different values of t. The closed circles refer to the 
X = case and the lines are a guide to the eye. 



both on the horizontal and vertical axis in the panels. In 
both panels, the closed circles correspond to the x=Q lat- 
tice case and the line is a guide to the eye. The results 
indicate a correlation between Tk and the f-c spin corre- 
lation function (the larger Tr-'s correspond to the more 
negative f-c values) as one would expect, since both pro- 
vide a measure of the Kondo effect. For t=0.05, most 
of the disordered cluster configurations are in the Kondo 
regime (Sg = 0), with larger Tk values; consequently, 
panel (a) has a larger number of singlet configurations. 
The introduction of MV impurities induces a distribution 
of Ti<-(i)'s, whose values are overall enhanced compared 
to those for the x=0 case, except for several configu- 
rations for t=0.05, in contrast with single-site theories 
for extended systems|^. It is interesting that P(Tk) for 
t—Q.Qb exhibits a bimodal behavior centered about 710 
and 570K, respectively: The higher Tr- 's originate from 
isolated Kondo atoms which have MV atoms as n.n. so 
that the local screening of the magnetic moment of the 
A atom is enhanced. 

The effect of alloying and A on the RKKY versus 
Kondo competition for a given x is seen in Fig. 4 
(left panel), where the configuration averaged local < 
Sf{i)S^{i) >^ and < Sf{i)Sf{i + 1) >^ correlation 
functions are plotted as a function of t. The solid curves 
denote the uniform x=0 case, where we drive the clus- 
ter from the RKKY to the Kondo regime as we decrease 
t. We find that the stronger the average Kondo correla- 
tions are the weaker the average RKKY interactions and 
vice versa. In the weak Kondo regime the configurations 
exhibit a wider distribution of RKKY interactions indi- 
cating that they are sensitive to the local environment. 
In contrast, in the strong Kondo regime, the Kondo (A) 
atoms become locked into local Kondo singlets and the 
n.n. RKKY interactions are insensitive to the local envi- 
ronment. Interestingly, both energy spacing and disorder 
lead to an overall enhancement of the RKKY interactions 
compared to the homogenous state. 

In the right panel of Fig. 4 we present the t versus x 
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FIG. 4: Left panel; Configuration-averaged local /-c (top) 
and n.n. /-/ spin correlations (bottom) for the A atoms as 
function of t. The solid line refers to the homogenous a; = 
case. Right panel: Zero-temperature t vs a; phase diagram 
for the nanocluster. Black (gray) circles denote the Kondo 
(RKKY) regime. The white circles and the dashed contour 
delimit the FM region. The horizontal stripes denote the non- 
magnetic (NM), mixed valence (MV) and local moment (LM) 
behavior of the B-atoms. 



phase diagram for the nanocluster at T = . We com- 
pare the < Sf{i)S^{i) >^ and < Sf(i)Sj(i + 1) >^ to 
assign a state of specific concentration to the Kondo or 
RKKY regimes (black and gray circles, respectively), in 
analogy with the a; = case (Table I) and with mean field 
treatments 24] . The horizontal gray stripes denote qual- 
itatively ranges of t where the B atoms exhibit NM, MV 
and LM behavior. An interesting feature of the phase 
diagram is the appearance of a large FM region (Sg ^ 0) 
enclosed by the dashed line. The RKKY region at large 
t and large x originates from the B atoms which become 
magnetic. For the non FM configurations and for a; < 5 
the Kondo (RKKY) correlations of the A atoms dominate 
at small (large) t, in analogy with the a; = case. On the 
other hand, for a; = 5 the local Kondo correlations of the 
single A atom at low t dominate over the /-/ correlations 
between the A-B and B-B pairs. For the uniform (a:=6) 
MV case we include only results in the large t regime, 
where the MV atoms acquire LM's which couple antifer- 
romagnetically. Overall, the RKKY interactions prevail 
for any concentration when t is comparable or larger than 
the hybridization V . 

In conclusion, we have presented exact diagonalization 
results for strongly correlated nanoclusters to study the 
effect of disorder and energy spacing on the interplay 
between the Kondo and RKKY interactions. Tuning A 
can drive the nanocluster from the Kondo to the RKKY 
regime, i.e. a tunable Doniach phase diagram in small 
clusters. The interplay of A and disorder produces a rich 
structure zero-temperature alloy phase diagram, where 
regions with prevailing Kondo or RKKY correlations al- 



ternate with domains of FM order. The distribution of 
local Tk and RKKY interactions depends strongly on 
the local environment and are overall enhanced by disor- 
der, in contrast to the hypothesis of single- impurity based 
"Kondo disorder" models for extended systems. The A 
versus disorder interplay may be relevant to experimen- 
tal realizations of small cluster with tunable magnetic 
properties. 
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